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O Abstract 

^^ We establish a connection between two previously unrelated topics: a particular dis- 

CN Crete version of conformal geometry for triangulated surfaces, and the geometry of ideal 

^^ polyhedra in hyperbolic three-space. Two triangulated surfaces are considered discretely 

Ci conformally equivalent if the edge lengths are related by scale factors associated with the 

^ vertices. This simple definition leads to a surprisingly rich theory featuring Mobius in- 

variance, the definition of discrete conformal maps as circumcircle preserving piecewise 

projective maps, and two variational principles. We show how literally the same theory 

can be reinterpreted to addresses the problem of constructing an ideal hyperbolic poly- 

1— I hedron with prescribed intrinsic metric. This synthesis enables us to derive a companion 

r^ theory of discrete conformal maps for hyperbolic triangulations. It also shows how the 

\^ definitions of discrete conformality considered here are closely related to the established 

(-H definition of discrete conformality in terms of circle packings. 
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Recall that tw^o Riemannian metrics g and g on a smooth manifold M are called conformally 
equivalent if 
ON g = e^"g (1.1) 

^ for a function u e C°°{M). In the discrete theory that w^e consider here, smooth manifolds are 

^ replaced with triangulated piecewise euclidean manifolds, and the discrete version of a confor- 

O rnal change of metric is to multiply all edge lengths with scale factors that are associated with 

'^ the vertices (Definition 2.1.1). Apparently, this idea to model conformal transformations in a 

• . discrete setting by attaching scale factors to the vertices appeared first in the four-dimensional 

.^ Lorentz-geometric context of the Regge calculus [32]. The Riemann-geometric version of this 

^ notion appeared in Luo's work on a discrete version of the Yamabe flow on surfaces [20]. 

He showed that this discrete Yamabe flow is the gradient flow of a locally convex function. 
Later, an explicit formula for this function was found (Ej.e.A defined in equation (3.3), with 
= 0), and this lead to an efficient numerical method to compute discrete conformal maps, 
suitable for applications in computer graphics [34]. (Much of the basic theory of conformal 
equivalence and conformal maps in Section 2 and the first variational principle in Section 3 
are already covered or at least touched upon in this earlier paper Since everything else in 
the present paper relies on this material, we deemed it necessary to review it here, including 
more mathematical detail than is appropriate in a graphics paper) The first variational prin- 
ciple described in Section 3.3 reduces the discrete conformal mapping problems (Section 3.1) 
to problems of convex optimization. Figures 1 and 2 show examples of discrete conformal 
maps that were obtained this way. Mobius transformations preserve the discrete conformal 





Figure 1. Discrete conformal map to a rectangle. 



class (Section 2.6), and this makes it possible to construct discrete conformal maps to regions 
bounded by circular polygons (Section 3.7), discrete analogs of the classical Riemann maps. 

The first variational principle involves a function of the (logarithmic) scale factors u. The 
second variational principle (Section 3.6) involves a function of the triangle angles. Both 
functions are related by a singular version of Legendre duality. The corresponding variational 
principles of the classical smooth theory are discussed in Section 3.8. 

There are clear signs in Sections 2 and 3 that indicate a connection w^ith hyperbolic geome- 
try: the appearance of Milnor's Lobachevsky function JI(x), the fact that the second variational 
principle is almost the same as Rivin's variational principle for ideal hyperbolic polyhedra with 
prescribed dihedral angles [29], and the definition of discrete conformal maps in terms of cir- 
cumcircle preserving piecew^ise projective functions (Section 2.7). This connection w^ith two- 
and three-dimensional hyperbolic geometry is the topic of Section 4. Reversing a construction 
of Penner [27] [7], we equip a triangulated piecewise euclidean surface with a canonical hy- 
perbolic metric with cusps. Discrete conformal maps are precisely the isometries with respect 
to this hyperbolic metric. The logarithmic edge lengths A (Section 2.1) and the length-cross- 
ratios that characterize a discrete conformal class (Section 2.3) are Penner coordinates and 
shear coordinates, respectively, of the corresponding hyperbolic surface (Section 4.3). The 
problem of flattening a triangulation discretely conformally is equivalent to constructing an 
ideal hyperbolic polyhedron with prescribed intrinsic metric (Section 4.4). With this inter- 
pretation of discrete conformality in terms of three-dimensional hyperbolic geometry, the two 
variational principles of Section 3 are seen to derive from Schlafli's differential volume formula 
and Milnor's equation for the volume of an ideal tetrahedron (Section 4.5). 

Once this connection between discrete conformality and hyperbolic polyhedra is estab- 
lished, it is straightforward to obtain a modified version of discrete conformality that per- 
tains to triangulations composed of hyperbolic triangles instead of euclidean ones (Section 5). 
This is the theory of discrete conformal uniformization of triangulated higher genus surfaces 
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Figure 2. Discrete Riemann maps. 



over the hyperbolic plane that has been applied for the hyperbolization of euclidean orna- 
ments [38]. (It is equally straightforward to obtain a corresponding theory for spherical tri- 
angulations, but the functions involved in the corresponding variational principles are not 
convex. We do not pursue this branch of the theory here.) 

The connection w^ith hyperbolic polyhedra entrains a connection between the discrete no- 
tion of conformality considered here and circle patterns, another discretization of the same 
concept. Thurston introduced patterns of circles as an elementary geometric visualization of 
hyperbolic polyhedra [36, Chapter 13]. He rediscovered Koebe's circle packing theorem [17] 
and showed that it followed from Andreev's work on hyperbolic polyhedra [1] [2]. Thurston's 
conjecture that circle packings could be used to approximate the classical Riemann map, which 
was later proved by Rodin and Sullivan [31], set off a flurry of research that lead to a full- 
fledged theory of discrete analytic functions and conformal maps based on packings and pat- 
terns of circles [35]. (The circle packing version of Luo's discrete Yamabe flow is the discrete 
Ricci flow of Chow and Luo [4].) The relationship between these two theories of discrete 
conformality is now clear: The circle packing theory deals with hyperbolic polyhedra with 
prescribed dihedral angles and the notion of discrete conformality considered here deals with 
hyperbolic polyhedra with prescribed metric. In Section 6.1 we discuss the relationship be- 
tween the variational principles for discrete conformal maps (Section 3) and two variational 
principles for circle patterns. One is due to Rivin [29] (see also the recent survey article by 
Futer and Gueritaud [11], which provides a wealth of material that is otherwise difficult to 
find), and the other is again related to it by the same sort of singular Legendre duality. 

In the last two sections, we discuss two extensions of discrete conformality that are moti- 
vated by the circle pattern connection discussed in Section 6.1. Instead of triangulations one 
can consider meshes composed of polygons that are inscribed in circles (Section 6.2), and we 
consider the problem to map multiply connected domains to domains bounded by polygons 
inscribed in circles, a discrete version of circle domains (Section 6.3). 

Two important questions are not addressed in this paper. The first is the question of con- 
vergence. Of course we do believe that (under not too restrictive assumptions that have yet to 
be worked out) discrete conformal maps approximate conformal maps if the triangulation is 
fine enough. Figure 2 clearly suggests that a version of the Rodin-Sullivan theorem [31] also 
holds in this case. But all this has yet to be proved. 

The other question concerns the solvability of the discrete conformal mapping problems of 
Section 3.1. A solution may not exist due to violated triangle inequalities. Fairly obvious nec- 
essary conditions and how they relate to properties of the function fj.e.A appearing in the first 
variational principle are discussed in Section 3.5. In the numerous numerical experiments that 
we have made, we have observed that a solution exists if the necessary conditions are satisfied, 
no triangles are almost degenerate to begin with, and the triangulation is not too coarse. But 
to find necessary and sufficient conditions for solvability seems to be an intractable problem 
in this setting. After all, this would amount to giving necessary and sufficient conditions for 
the existence of a (not necessarily convex) ideal hyperbolic polyhedron with prescribed in- 
trinsic metric and prescribed combinatorial type. The way out is to restrict oneself to convex 
polyhedra while widening the concept of discrete conformal map to allow for combinatorial 
changes (Section 4.1). Rivin proved that any hyperbolic metric with cusps on the sphere is 
realized by a unique ideal polyhedron [30]. This translates into an existence statement for 
discrete conformal maps. We will discuss the issue of solvability in greater detail in a future 
paper [33], in which we will present a constructive proof for Rivin's realization theorem that 
is based on the methods presented here. (The approach is very similar to the recent construe- 



tive proof of Alexandrov's classical polyhedral realization theorem [3].) Rivin considered the 
non-constructiveness of his proof its greatest drawback [30]. Although the theorem has in the 
meantime been greatly generalized [8] [9], all proofs are non-constructive. 

2 Discrete conformal equivalence and discrete conformal maps 

2.1 Discrete conformal equivalence 

A surface is a connected 2-dimensional manifold, possibly with boundary. A surface triangu- 
lation, or triangulation for short, is a surface that is a CW complex whose faces (2-cells) are 
triangles which are glued edge-to-edge. We will denote the sets of vertices (0-cells), edges (1- 
cells), and faces of a triangulation T by Vj, Ej, and Tj, and we will often drop the subscript 
T if the triangulation is clear from the context. We will also write Aj for the set of triangle 
angles, where angles means corners, or triangle-vertex incidences, not angle measures. 

A eucUdean surface triangulation, or euclidean triangulation for short, is a surface trian- 
gulation equipped with a metric so that T \ Vj is locally isometric to the euclidean plane, 
or half-plane if there is boundary, and the edges are geodesic segments. In other words, a 
euclidean surface triangulation is a surface consisting of euclidean triangles that are glued 
edge-to-edge. At the vertices, the metric may have cone-like singularities. 

A euclidean triangulation is uniquely determined by a triangulation T and a function i : 
Ej — > M>o assigning a length to every edge in such a way that the triangle inequalities are 
satisfied for every triangle in Tj. We call such a positive function i on the edges that satisfies all 
triangle inequalities a discrete metric on T, and we denote the resulting euclidean triangulation 
by(T,£). 

In this paper, we will assume for simplicity that the triangulations are simplicial complexes. 
This means that a triangle may not be glued to itself at a vertex or along an edge, and the 
intersection of two triangles is either empty or it consists of one vertex or one edge. This 
restrictions to simplicial complexes allows us to use simple notation: we will denote by ij the 
edge with vertices i and j, by ijk the triangle with vertices i, j, and k, and by j^. the corner 
at vertex i in triangle ijk. li f,g,h, and (j) are functions on V, E, T, and A, respectively, we 
will vwite /;, gij, hij]^, and (^'j, for /(i), g(ij), h(ijk), and 0(jfc). But while this restriction to 
simplicial complexes is notationally very convenient, it is a priori uncalled for. There are a few 
exceptions, like Section 3.7 on mapping to the sphere and disk, but in general the domain of 
validity of the theory presented here extends beyond the simplicial case. (This will be essential 
in the sequel [33].) 

2.1.1 Definition (Luo [20]). Two combinatorially equivalent euclidean triangulations, (T, £) 
and (T, I), are discretely conformally equivalent if the discrete metrics I and £ are related by 

£.. ^gi(u,+Uj)^.. (2.1) 

for some function u : V — > M. This defines an equivalence relation on the set of discrete metrics 
on T, of which an equivalence classes is called a discrete conformal class of discrete metrics, or 
a discrete conformal structure on T. 

Instead of the edge lengths £ we will often use the logarithmic lengths 

A = 21og£. (2.2) 




Figure 3. The length-cross-ratio on edge ij. The lengths of the bold solid and bold dashed edges 
appear in the numerator and denominator, respectively. 



(The reason for the factor of 2 will become apparent in Section 4.) In terms of these logarithmic 
lengths, relation (2.1) between I and £ becomes linear: 



^ij = '^ij + "; +Uj. 



(2.3) 



2.2 The two most simple cases 



(1) If the triangulation T consists of a single triangle ijk, then any two euclidean triangula- 
tions (T, £) and (T, £) are discretely conformally equivalent, because the three equations 



always have a unique solution for Uj, Uj and Uj,: 
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(2) Now let T be the triangulation consisting of two triangles ijk and ilj glued along edge 
ij as shovwi in Figure 3, and let (. and I be two discrete metrics on T. What is the condition 
for (T, £) and (T, £) to be discretely conformally equivalent? For each triangle considered 
separately, the corresponding equations (2.1) determine unique solutions for the values of 
u on its vertices. For each of the common vertices i and j one obtains two values and the 
necessary and sufficient condition for discrete conformal equivalence is that they are equal, 
which is equivalent to the condition 
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2.3 Length-cross-ratios 

The simple reasoning of Section 2.2 extends to the general case: Let T be any triangulation, 
and let I and i be two discrete metrics on T. For each triangle ijk of T, considered separately, 
equations (2.1) determine unique values for u. Thus, for each vertex i e V, one obtains one 
value for Uj per adjacent triangle. These values are in general different. They agree for each 
vertex if and only if the discrete metrics £ and i are discretely conformally equivalent. Since 



the vertex links are connected, it suffices to consider values obtained from adjacent triangles. 
This leads to Proposition 2.3.2 belowr, w^here the condition for discrete conformal equivalence 
is given in terms of the so-called length-cross-ratios: 

2.3.1 Definition. For each interior edge ij betw^een triangles ijk and ilj as in Figure 3, define 
the length-cross-ratio induced by i to be 

This definition implicitly assumes that an orientation of the triangulated surface has been 
chosen. The other choice of orientation leads to reciprocal values for the length-cross-ratios. 
(For non-orientable surfaces, the length-cross-ratios are w^ell defined on the interior edges of 
the oriented double cover.) 

If the quadrilateral iljk is embedded in C, then the length-cross-ratio Ictjj is just the abso- 
lute value of the complex cross ratio of the vertex positions Zi,Zi,Zj,z^, 

, . (^1 ~ ^2)(Z3 - Z4) 

cr(zi,Z2,Z3>Z4) = J- —rr rr • 

(Z2-Z3XZ4-Z1J 

Discretely conformally equivalent metrics I, l induce the same length-cross-ratios, because 
the scale factors e"'^ cancel. By the reasoning above, the converse is also true. 

2.3.2 Proposition. Two eucUdean triangulations (T, £) and (T, £) are discretely conformally 
equivalent if and only if for each interior edge ij e Ej, the induced length-cross-ratios are equal: 

2.4 The product of length-cross-ratios around a vertex 

Let us denote the sets of interior edges and interior vertices by £j„f and Vj„(., respectively. Which 
functions £j„f — > M>o can arise as length-cross-ratios? A necessary condition is that the product 
of length-cross-ratios on the edges around an interior vertex is 1, because all lengths I cancel. 

For all i e V^^^ : ]~[ Icr;^ = 1 . (2.5) 

j:ijeE 

If we ignore the triangle inequalities, this condition is also sufficient: 

2.4.1 Proposition. Let Icr : Ej^f — > M>o be any positive function on the set of interior edges. There 
exists a positive function £ : £ — > M>o on the set of edges satisfying (2.4) for every interior edge ij, 
if and only if condition (2.5) holds. 

Proof. It remains to show^ that if Icr e (M>o)^'"' satisfies condition (2.5), then the system of 
equations (2.4) has a solution. In fact, w^e w^ill explicitly construct such a solution. To this 
end, w^e introduce auxiliary parameters c, which are defined on the set of angles A of the 
triangulation: Given I e (M>o)^, define c e (M>o)^ by 

see Figure 4. In terms of these parameters, the length-cross-ratios induced by £ are 




I 
Figure 4. The parameters c^.^, defined on the set of triangle angles A 

d 
lcr;, = ^, (2.7) 

where I, j, k occur in the link of i in this cyclic order, as in Figure 3. (For a geometric 
interpretation of the parameters c'^^ in terms of hyperbolic geometry, see Section 4.2.) 

Now suppose Icr e (M>o)^'"' satisfies condition (2.5). Then it is easy to find a solution c e 
(M>o)^ of equations (2.7), because each equation involves only two values of c on consecutive 
angles at the same vertex. So one can freely choose one c-value per vertex and successively 
calculate the values on neighboring angles around the same vertex by multiplying (or dividing) 
with the values of Icr on the edges in between. 

Next, solve equations (2.6) for I, where c is the solution to equations (2.7) just constructed. 
This is also easy: The length of an edge ij is determined by the values of c on the two adjacent 
angles on either side, 

(Check that the two c-values on the other side give the same value.) Thus we have constructed 
a function I e (M^g)^ satisfying equations (2.4) for the given function Icr e (M^q)^'"'^- □ 

2.5 Some spaces and dimension counts 

The set of discrete metrics on a triangulation T is the open subset of M^ containing all £ e M^ 
that satisfy £ > and the strict triangle inequalities. Since it is not empty (always containing 
the constant function 1^) and the inequalities are linear, it is the interior of an |£| -dimensional 
convex polytope. The discrete conformal structures on T are submanifolds which form a fi- 
bration of the polytope of discrete metrics. If £ is a discrete metric, its conformal class is 
parametrized by u e M^ via equation (2.1). The space of discrete conformal classes has every 
right to be called the discrete Teichmilller space of the triangulation T. It is parameterized by 
the length-cross-ratios. But these are excessive parameters, for they are constrained by the 
product conditions (2.5), see Section 2.4. 

The logarithmic lengths A defined by equation (2.2) are coordinates on the space of dis- 
crete metrics. The space of discrete metrics corresponds to the parameter domain A c M^ of 
all A e M^ satisfying 

-e^y/^ + e^i'/2 + e^^^l'^ > 

gAy/2 _ gA,t/2 _^ gA;,,/2 ^ Q (-28) 
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for all triangles ijk e T. While the domain A is not a polytope any more, not even a convex 
set, the fibration by conformal classes has become simpler. For A e M^ let U^ be the affine 
subspace of all A parameterized by u e M^ as in equation (2.3). If A e A, its discrete conformal 
class is U^nA. If A ^ A, the intersection may be empty; if it is not empty, it is also a discrete 
conformal class. The fibration of the space of discrete metrics by conformal classes corresponds 
to a fibration of the domain A by the family {[/aIasm^ of parallel affine subspaces, and the 
discrete Teichmiiller space corresponds to an open set in the factor space M^/t/g- Thus, we 
have: 

2.5.1 Proposition. The dimension of each discrete conformal class is \Vj\, and the dimension of 
the discrete Teichmiiller space is \Ej\ — \Vj\. 

Example. Consider the case where T is a triangulation of the closed orientable surface of 
genus g. Using2-2g = |r|-|£| + |V| and 3|r| = 2|£|, one obtains |£| - |y| = 6g -6 + 2|V|. 
This is also the dimension of 5^^ |y|, the Teichmiiller space of a genus g Riemann surfaces with 
|V| punctures. It is no coincidence that the dimensions agree. The discrete conformal classes 
actually correspond to points in the Teichmiiller space ^g^\v\ (see Section 4.1). 

The subspaces U^ were described above in parameter form. An equivalent description in 
equation form is also available. If Icr : M^'"' — > M>o is the length-cross-ratio function induced 
by £ = e'*-/^, then U;>^ is the space of all A e M^ satisfying 

^ii - hj + ^jk - hi = 21oglcr;j- (2.9) 

for all ij e £j„f, where k and I are the opposite vertices as in Figure 3. The logarithmic length- 
cross-ratios log Icr are coordinates for the discrete Teichmiiller space; but they are excessive 
coordinates because they are constrained by the condition (2.5) that their sum around each 
interior vertex is zero. From this one obtains again |£| — |y| (#variables — #constraints) for 
the dimension of the discrete Teichmiiller space. 

2.6 Mobius in variance of discrete conformal structures 

The group of Mobius transformations of M" = M" U |oo} is the group generated by inversions 
in spheres. (Planes are considered spheres through oo.) The group of Mobius transformations 
is also generated by the similarity transformations (which fix oo), and inversion in the unit 
sphere. Mobius transformations are conformal, and a famous theorem of Liouville says that 
for n > 2, any conformal map of a domain [/ c M" is the restriction of a Mobius transformation. 
Let T be a triangulation and let || • || denote the euclidean norm on M", n > 2. Suppose 
V : Vj — > M" maps the vertices of each triangle to three affinely independent points. Then v 
induces a discrete metric ^ij = ||Vj — Vj||. Two maps v, i/ : V — > M" c M" are related by a Mobius 
transformation if there is a Mobius transformation T such that v = T ov. 

2.6.1 Proposition. If two maps Vj — > M" are related by a Mobius transformation, then the 
induced discrete metrics are discretely conformally equivalent. 

Proof The claim is obvious if the relating Mobius transformation is a similarity transformation. 
For inversion in the unit sphere, x >-> ^— ^ x, it follows from the identity 
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Figure 5. A coarsely triangulated domain in the plane (middle) is mapped to a rectangle (right) 
by a discrete conformal map (see Definition 2.7. 1). Instead of using circumcircle preserving piece- 
wise projective interpolation, one can also interpolate linearly in each triangle. The result (left) 
looks noticeably "less smooth". 

2.6.2 Remark. For n = 2 there is an obvious alternative argument involving the complex cross 
ratio. One can extend this argument to n > 2. The only difficulty is to define a complex cross- 
ratio for four points in M" if n > 2, such that it is invariant under Mobius transformations. Such 
a cross-ratio can be defined up to complex conjugation by identifying a 2-sphere through the 
four points conformally with the extended complex plane C. This involves several choices: a 
choice of 2-sphere if the four points are cocircular, a choice of orientation of the 2-sphere, and 
choice of orientation preserving conformal map to C. Only the choice of orientation makes a 
difference, the two choices leading to conjugate values for the cross-ratio. The length-cross- 
ratio is the absolute value of this complex cross-ratio, so the ambiguity with respect to complex 
conjugation does not matter. 

2.7 Discrete conformal maps 

So far we have only talked about discrete conformal equivalence. This section deals with the 
matching notion of discrete conformal maps (see Figure 5) . 

For any two euclidean triangles (with labeled vertices to indicate which vertices should be 
mapped to which), there is a unique projective map that maps one triangle onto the other and 
the circumcircle of one onto the circumcircle of the other (see Lemma 2.7.3 below). Let us call 
this map the circumcircle preserving projective map between the two triangles. 

2.7.1 Definition. A discrete conformal map from one euclidean triangulation (T, £) to a com- 
binatorially equivalent euclidean triangulation (T, £) is a homomorphism whose restriction 
to every triangle is the circumcircle preserving projective map onto the corresponding image 
triangle. 

Consider two combinatorially equivalent euclidean triangulations, (T, £) and (T, £). For 
each individual triangle of (T, £), there is a circumcircle preserving projective map to the cor- 
responding triangle of (T, £). But these maps do in general not fit together continuously across 
edges. However, they do fit together, forming a discrete conformal map, precisely if the eu- 
clidean triangulations are discretely conformally equivalent: 
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2.7.2 Theorem. The following two statements are equivalent: 
(i) (T, £) and (T, £) are discretely conformally equivalent. 

(ii) There exists a discrete conformal map (T, £) — > (T, £). 

The rest of this section is concerned with the proof of Theorem 2.7.2. It follows easily from 
Lemma 2.7.3 below, which provides an analytic description of the circumcircle preserving 
projective map between two individual triangles. 

Consider two triangles A and A in the euclidean plane, and let (x;,/;) and (x;,/;), i e 
{1,2,3}, be the coordinates of their vertices in a Cartesian coordinate system. Let i^ and i^j 
be the side lengths, 

and similarly for iij. Consider the euclidean plane as embedded in the projective plane MP^ 
and let Vj = (x^,y^,l) and Vj = (Xj,yj, 1) be the homogeneous coordinate vectors of the ver- 
tices, normalized so that the last coordinate is 1 . Then the projective maps / : MP^ — > RP^ that 
map A to A correspond via /([v]) = [f (v)] to the linear maps f : M"' — > M"' of homogeneous 
coordinates that satisfy 

f(v;) = M:Vi (2.10) 

for some "weights" jU; e M \ {0}. 

2.7.3 Lemma. The projective map f : [v] >-> [f (v)] maps the circumcircle of A to the circumcir- 
cle of A if and only if 

(/ii,M2,/i3) = M(e""\e""^e""^) (2.11) 

where Ui,U2>"3 ^^^ the logarithmic scale factors satisfying the three equations (2.1) for a single 
triangle and ;U e M \ {0} is an arbitrary factor 

Proof of Lemma 2.7.3. The circumcircle of A is 

{[v]eMp2|q(v) = 0}, 

where q is the quadratic form 

q(x,y,z) = x^ -l-y^ -|-2axz + 2byz + cz^ 

with a, b, c e M uniquely determined by the condition that 

q(vi) = q(v2) = q(v3) = 0. (2.12) 

In the same fashion, let the quadratic form describing the circumcircle of A be 

q(x,y,z) = x^ + y^ + 2dxz + 2hyz + cz^ 

so that 

q(vi) = q(v2) = q(v3) = 0. (2.13) 

We will also denote by q and q the corresponding symmetric bilinear forms: 

q(v) = q(v,v) and q(v) = q(v,v). 
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The projective map / maps circumcircle to circumcircle if and only if q and the pull-back 
F*q are linearly dependent. That is, if and only if 

/i^ q(v, w) = q(F(v), F(w)) 

for all u, w e M'^ for some /i e M. Since v^, V2, V3 is a basis of M"^ and because of equations (2.12) 
and (2.13), this is the case if and only if 

M^q(Vi, Vj) = /i;/ljq(V;,Vj) (2.14) 

for i,j e {1,2,3}, i j^ j. Now note that 

^l = q(Vi - Vj, V; - Vj) = -2q(v;, v^, 
and similarly £^. = —2q(vi,Vj). So condition (2.14) is equivalent to 

Solve equations (2.1) for U; to obtain (2.11). This completes the proof of Lemma 2.7.3. D 

To prove Theorem 2.7.2, consider two euclidean triangulations (T, £) and (T, £), and a pair 
of adjacent triangles ijk and jil of T. Embed the corresponding euclidean triangles of (T, £) 
simultaneously isometrically in the euclidean plane, and do the same for the corresponding two 
euclidean triangles of (T, £). Lemma 2.7.3 tells us what the circumcircle preserving projective 
maps are, and we might as well choose ju = 1 in both cases. These two maps fit together 
continuously along edge ij if and only if the values of /i, = e~"' and [jLj = e~"j from one 
triangle are proportional to those of the other triangle. Since the value of e'^"'+"i^/^ = lij/lij 
is the same for both triangles, this is the case if and only if the values of ji^ and jUj, hence 
also those of Ui and Uj, coincide for both triangles. This holds for all interior edges if and 
only if (T, £) and (T, £) are discretely conformally equivalent. This completes the proof of 
Theorem 2.7.2. 

3 Two variational principles 

3.1 Discrete conformal mapping problems 

Consider the following type of discrete conformal mapping problem: 

3.1.1 Problem. Given 

• a surface triangulation T, 

• a discrete conformal class "iS of discrete metrics on T, 

• a desired angle sum 0, for each vertex i e Vj, 

find a discrete metric £ in the conformal class "iS such that the euclidean triangulation (T, £) 
has angle sum 0j around each vertex i e Vj. 

In particular, if the given desired angle sum 0, equals 2n for every interior vertex i, then 
Problem 3.1.1 asks for aflat euclidean triangulation in the given conformal class which has 
prescribed angles at the boundary. A flat and simply connected euclidean triangulation can 
be developed in the plane by laying out one triangle after the other. Thus, Problem 3.1.1 
comprises as a special case the following problem. 

12 



3.1.2 Problem. Given 

• A euclidean triangulation (T, £) that is topologically a disc and 

• a desired angle sum 0, for each boundary vertex i, 

find a discretely conformally equivalent planar triangulation (T, £) w^ith the given angle sums 
at the boundary. (Self-overlap is allow^ed.) 

We also consider a more general type of problem than Problem 3.1.1. Suppose the discrete 
conformal class Sg' is given in the form of a representative metric £ e (M>o)^. For some vertices 
i we may prescribe the (logarithmic) scale factor U; instead of the angle sum 0;: 

3.1.3 Problem. Given 

• a triangulation T, 

• a function I e (M>o)^ representing a conformal class, 

• a partition V = Vg U V^ of the vertex set, 

• a prescribed logarithmic scale factor U; e M for each vertex i e Vq 

• a prescribed angle sum 0; for each vertex i e V^ 

find logarithmic scale factors U; e M for the remaining vertices i e V^ so that I determined by 
equations (2.1) is a discrete metric and (T, £) has the prescribed angle sum 0j around each 
vertex i e V^ . 

For Vg = 0, Vj = V, this is just Problem 3.1.1. (If the conformal class S^ is given in the 
form of length-cross-ratios Icr e (M>o)^'", one can obtain a representative I e (M>o)^ using 
the method described in the constructive proof of Proposition 2.4.1.) Note that any instance 
of Problem 3.1.3 can be reduced to the special case w^here U; = is prescribed for i e Vg- 
Simply apply first a discrete conformal change of metric (2.1) with the arbitrary prescribed Uj 
for i e Vg. 

In Section 3.7 w^e w^ill see how^ one can also find discrete conformal maps to the sphere 
and to the disk (that is, solve the discrete Riemann mapping problem) if one can solve Prob- 
lem 3.1.3. 

3.1.4 Remark. An important special case of Problem 3.1.3 is the following: prescribe the 
angle sum 0; = 271 for interior vertices, and U; on the boundary. This is analogous to the 
following boundary value problem of the smooth theory: Given a smooth 2-manifold with 
boundary M equipped with a Riemannian metric g, find a conformally equivalent flat Rie- 
mannian metric e^"g with prescribed uI^m- Suppose we measure the relative distortion of a 
conformally equivalent Riemannian metric by the Dirichlet energy of u, D(u) = - J du A *du. 
Then the conformally equivalent flat Riemannian metrics with least distortion are those with 
^\dM — const. Thus, up to scale there is a unique least distortion solution obtained by setting 
"Ism = 0- [34, Appendix E] 

3.1.5 Theorem. If Problem 3.1.3 (with Vg ^ @) has a solution, it is unique. If Problem 3.1.1 has 
a solution, it is unique up to scale. In any case, if a solution exists, it can be found by minimizing 
a convex function. 

Proof. This follows from Propositions 3.3.3, 3.3.4, and 3.3.5. D 

3.2 Analj^ic description 

Analytically, these discrete conformal mapping problems amount to solving systems of nonlin- 
ear equations. Suppose £ e S^ is a representative of the given discrete conformal class, and 
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consider £ as function of u defined by equations (2.1). Then solving the general Problem 3.1.3 
amounts to solving for the unknown U; (i e V^) a system of angle-sum equations 

S a;, = e,, (3.1) 

}k:ijkeT 

with one equation for each vertex t e V^. Here, a'j^. is the angle at i in triangle ijk of (T, £). 

The angles a are nonlinear functions of the new lengths i. They can be obtained by invoking, 
for example, the cosine rule or the half-angle formula 




+ ^jk + hiXhj + ijk-hi) 



2 ; V <i^ij-^jk+ki)(^ij+^jk+hi) 

(Tilde marks over a and i have been omitted in this equation to avoid visual clutter) 

3.3 The first variational principle 

The system of nonlinear equations described in the previous section turns out to be varia- 
tional. Solutions of the conformal mapping problems correspond to the critical points of the 
function fj.e.A defined as follows. The precise statement of this first variational principle is 
Proposition 3.3.3. 

Let T be a surface triangulation, 9 e M^, and A e M^. For now (we will later extend the 
domain of definition to M^) define the real valued function Ej.o.aC") on the open subset of 
M^ containing all u such that £ determined by (2.1) is a discrete metric (that is, satisfies the 
triangle inequalities) : 



Et,0,x{u) = Y, («o-^0- + ^ikhk + 4i~^ki + 2JI(af.) + 2JI(a;.,) + 2JI(ai;) 

ijkeT 



The first sum is taken over all triangles, a'j, denotes the angle at vertex i in triangle ijk with 

side lengths I = e^^^, X is the function of u defined by equations (2.3), and JI(x) is Milnor's 
Lobachevsky function, 

JI(x) = - log|2sin(t)|dt. (3.4) 

Jo 
(Figure 6 shows a graph of this function.) The second sum is taken over all vertices. It is linear 
in u. 

3.3.1 Remark. The notation JI(x), using a letter from the Cyrillic alphabet, and the name 
"Lobachevsky function" are due to Milnor [22] [23]. Lobachevsky used a slightly different 
function which is also known as the Lobachevsky function and often denoted L(x). To distin- 
guish these two functions, we call JI(x) Milnor's Lobachevsky function. It is almost the same 
as Clausen's integral [19], ClaCx) = |jl(2x). 

3.3.2 Proposition (First derivative). The partial derivative of Ejqx with respect to u^ is 

' jk:i]keT 

where the sum is taken over all angles around vertex i. 
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Figure 6. Graph of Milnor's Lobachevsky function, y = JI(x). The function is n-periodic, odd, 
and smooth except at x ^ nZ, where its tangents are vertical. 



Proof. This follows from equation (3.7) and Proposition 3.4.1. 



D 



3.3.3 Proposition (First variational principle). Let "iS be a discrete conformal class with repre- 
sentative metric I = e^^^, and let I = e^^^ where A is the function ofu defined by equations 2.3. 
Then 

• l solves Problem 3.1.1 if and only if u is a critical point of Ejq^^ 

• l solves Problem 3.1.3 if and only if u is a critical point of Ejqx with fixed u^for i e Vq. (In 
this case, the values of Q for i e Vq cl^^ irrelevant.) 

rs 

Proof. This follows immediately from Proposition 3.3.2 because ^ ^T,e,A = is equivalent to 
the angle sum condition (3.1). D 



3.3.4 Proposition (Local convexity). The function £t,0,a ^•^ locally convex, that is, its second 
derivative ^ "^^'^ du; duj is a positive semidefinite quadratic form. The kernel is 1-dimensional 
and consists of the of the constants in M^. 

Proof. This follows from equation (3.7) and Proposition 3.4.3. D 

3.3.5 Proposition (Extension). The function £7,0,1 can be extended to a convex continuously 
differentiable function on M^. 



Proof. This follows from equation (3.7) and Proposition 3.4.4. 



D 



In fact, one has an explicit formula for the second derivative of Ej.e.A- This is helpful from 
the practical point of view, because it allows one to use more powerful algorithms to minimize 
EjQx and thus solve the discrete conformal mapping problems. It is also interesting from 
the theoretical point of view, because the second derivative of £T,e,A ^t u is the well known 
finite-element approximation of the Dirichlet energy (the cotan-formula) for a triangulation 
with edge lengths I [6] [28]: 

3.3.6 Proposition (Second derivative). The second derivative of Ejqx at " is 

S a a'"' ^"i rf"i =0X1 w; •(u)(du; - dUjf, 
. . „ dUjdu; ■' 2 ^—L ■' ■' 
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Figure 7. Triangle with sides a = e^, b = e^, c = e^ and angles a, 13, y. The radius of the 
circumcircle is R= —4- 



2 sin a 2sin/3 2 sin 7 

where WiAu) = -(cotaf. + cota|.) if ij is an interior edge with opposite vertices k and I and 

J ^ '■J '■J 

Wij{u) = -cota*^. if ij is a boundary edge with opposite vertex k. This assumes all triangle 
inequalities are satisfied. If triangle inequalities are violated, the cotangent terms for the corre- 
sponding triangles have to be replaced with 0. 

Proof. This follows from equation (3.7) and Proposition 3.4.2. D 

3.4 A peculiar triangle function 

Consider the function 

/(x, y,z) = ax + Py + rz + JI(a) + JI(^) + JI(r), (3.5) 

where a, p, and y are the angles in a euclidean triangle with sides a = e^, b = e^ , and c = e^ 
as shown in Figure 7. Such a triangle exists if and only if the triangle inequalities are satisfied. 
So / is (for now) only defined on the set 

j^ = {(x,y,z)eM^| -e^' + e^ + e^ >0, 

e^ - e^ + e" > 0, (3.6) 

e^ + e^-e">0}. 

The function /(x, y, z) is the fundamental building block of Ejq^(u) since 

ET,e,xM= Xi [2fi^,^,'f)-^iXij,Xj,,X,i))+J]QiUi, (3.7) 

ijkeT ieV 

and Propositions 3.3.2, 3.3.4, 3.3.5, and 3.3.6 follow from corresponding statements regarding 

f(.x,y,z). 

3.4.1 Proposition (First derivative). The partial derivatives off are 

df df df 

ox oy oz 

Proof. Using JI'(x) = — log |2 sin(x)| we obtain from (3.5) that 

df ^ Sa ^ 5/3 ^ S7 

-- = a+ [x - log(2sma)J V [y - log(2sm/3)J V [z - log(2smy-)J -- . 

ox ox ox ox 
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Since 

X — log(2 sin a) = y — log(2 sin ^) = z — log(2 sin 7) = logR, 

where R is the radius of the circumcircle, and since 

da dB dr 

— + -^ + — ^ = 
dx dx dx 



3f 



a. 



(because a + ^ + 7 = 71), this implies 

3.4.2 Proposition (Second derivative). The second derivative of f is 



U 



[dy] 



f £^ 3^ fLL\ 

dx'^ Sx3y Sx3z 

d^f d^f a^f 

dydx 5y2 dydz 

V dzdx dzdy dz^ J 



\dy\ = cot a (dy - dz)'^ + cot P [dz - dx)^ 

^ dz -^ 

+ cot 7 (dx-dy)^ (3.8) 



Proof. By Proposition 3.4.1, 

(") 



I 3x2 3x3y 3x3z ' 

a^/ a\f_ 8^f 

3ydx 3y2 3ydz 

■ 3^ a^ 3jj_ . 

V dzdx dzdy d-^ J 



= {da dl5 dy), 



so the left-hand side of equation (3.8) equals 

dadx-l-d/3dy + dYdz. 

We w^ill show^ that 

da = cot7(dx — dy) + cot 13 (dx — dz). 

This and the analogous equations for d/3 and dy imply 

dadx-|-d/3dy-|-d7dz = cota(dy - dz)^ + cot^(dz - dx)^ + cot7(dx - dy)^ 

and hence equation (3.8). 

To derive equation (3.9), differentiate the cosine rule 

2bc cos a= b^ + c^ - a^ 

to get 

— 2i)c sin a da -|-2bc cos a (dy + dz) = 2b^ dy + 2c^dz — 2a^dx. 

Apply the cosine rule three more times to get 

2bc sin ada = (b^ - c^ + a^)(dx - dy) + (-b^ + c^ + a^)(dx - dz) 
= 2abcos7(dx — dy) + 2ac cos (3 (dx — dz) . 

Divide through by 2bc sin a and apply the sine rule to obtain equation (3.9). 



(3.9) 



D 



3.4.3 Proposition (Local convexity). The function f is locally convex, that is, the second deriva- 
tive (3.8) is a positive semidefinite quadratic form. Its kernel is one-dimensional and spanned by 

(l,l,l)eM3. 
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Figure 8. Left: Contour plot of (x,y) >-> /(x,y, 0). The intersection of the domain j4 (see 
equation (3.6)j with the xy-plane is shaded. Right: Graph of the same function, also showing 
contour lines. 



Proof. Writing [dy — dz) as ((dy — dx) — (dx — dz)) we obtain 

cot a {dy - dz)'^ + cot 13 {dz - dx)^ + cot y (dx - dy)^ 

= (cota + cot^)(dx - dz)^ + (cota + cot'y-)(dx -dy)^ 

— 2 cot a(dx — dy)(dx — dz) 

Thus, in terms of (dx — dz) and (dx — dy), the matrix of this quadratic form is 



M 

Using a + p + Y = n, we obtain 

1 



M 



cot a + cot p — cot a 
— cot a cot a + cot y 



smy 



— cos a sin ^ sin y 
sin a sin /3 sin 7 I -cos a sin /3 sin y sin^ 



and detM = 1. Since M^ > and detM > 0, M is positive definite. The claim about the 
second derivative of/ follows. D 

3.4.4 Proposition (Extension). Extend the definition of f from ^ to R as follows. Define 
/(x,y, z) by equation (3.5) for all (x,y,z) e M.^, where for (x,y,z) ^ j^ the angles a, /3, and y 
are defined to be n for the angle opposite the side that is too long and for the other two. The so 
extended function / : M'^ — > M is continuously differentiable and convex. 

Figure 8 shows contour lines of the extended function /(x, y,z) in the plane z = and its 
graph. 

Proof. The so-defined functions a, (3, y are continuous on M"^. This impHes the continuity of 
/ and, together with Proposition 3.4.1, the continuity of its first derivative. Since / is locally 
convex in j?/ (Proposition 3.4.3) and linear outside, it is convex. D 

3.4.5 Remark (Amoebas and Ronkin functions). In fact, js/ is an amoeba, the extended / is 
a Ronkin function, and the convexity of / that we have proved by elementary means follows 
also from a general theorem of Passare and Rullgard [26], which says that a Ronkin function 
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is convex. Amoebas were introduced by Gelfand, Kapranov and Zelevinsky [12]. The amoeba 
j4p of a complex polynomial p(zi, . . . ,Zn) with n indeterminates is defined as the domain in M" 
that is the image of the set of zeros of p under the map (z^,. ..,2^) >-> (log|zi|,...,log|z„|). So 
the domain j4 defined by equation (3.6) is the amoeba of the linear polynomial 2^+22 + 23. 
The Ronkin function of a polynomial p is defined as the function iV^ : M" — > M, 

if I I dzj dz^ 

JVp(xi,...,xJ=— — log|p(zi,...,zJ|— A...A— , 

where S^(r) is the circle in C around with radius r. As it turns out, 

We will not spoil the reader's fun by presenting a proof here. The same Ronkin function also 
appears in the work of Kenyon, Okounkov, and Sheffield on the dimer model [15] [16] as 
the Legendre dual of a "surface tension" in a variational principle governing the limit shape 
of random surfaces. (See in particular Kenyon's survey article on dimers [14] and Mikhalkin's 
survey article on amoebas [21].) Whether or how this is related to the variational principles 
discussed in this paper is unclear. 

3.5 Necessary conditions for the existence of a solution 

In this section, we will discuss some rather obvious and rather mild necessary conditions for 
the solvability of the discrete mapping problems and how they relate to the behavior of the 
function Ej.b.aC")- In short, the conditions are necessary for the problems to have a solution 
and sufficient to ensure that Ej.g.aC") behaves "sanely", so that the following solvability al- 
ternative (see the corollary to Proposition 3.5.3) holds: Provided that we are able to find a 
minimizer of a convex function if it exists, then the variational principle allows us to either 
solve a discrete conformal mapping problem or to ascertain that it is not solvable. 

Condition A. Siev ®! = ^ I ^ I 

If Problem 3.1.1 has a solution then clearly Condition A is satisfied (because the sum of angle 
sums around vertices equals the sum of angle sums in triangles). This is actually a discrete 
version of the Gauss-Bonnet formula. If we set K^ = 271—0; for interior vertices and Kj = ti— 0; 
for boundary vertices then Condition A is equivalent to 

Y^^i+Yj '^; = 2<|r|-|£| + |V|)). 
3.5.1 Proposition. The function £t,0,a(") ^ scale-invariant, that is, 

if and only if Condition A is satisfied. 

Proof Because /(x + h,y + h,z + h) = f(x,y,z) + nh , 
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The following stronger Conditions B and C are also obviously necessary for the existence of 
a solution of Problem 3.1.1. Moreover, if a solution to the general Problem 3.1.3 exists (w^here 
0j is given only for i e Vi), then positive ©-values can be assigned also to the vertices in Vq so 
that Conditions B and C are satisfied. 

Condition B. There exists a system of angles a > 0, such that 

a'j, + a^. + a. . = 71 for all ijk^T, 

and 

^ ajfc = 0i for all i e V. 

jk:ijkeT 

Condition C. If T^ is any subset of T and Vj c y is the set of all vertices of the triangles in T^, 
that is, 

Vi= U {i,j,k}, 

ijkeT 

then 

n\T\T^\> Y, ®i' 

ieV\V-^ 

where equality holds if and only if T^ = or T^ = T. 

3.5.2 Proposition. Conditions B and C are equivalent. 

Proof. The implication 'Condition B => Condition C is easy to see. Regarding the converse 
implication, Colin de Verdiere proves a similar statement using the feasible flow theorem [5, 
Section 7]. It is straightforward to adapt his proof for Proposition 3.5.2. D 

3.5.3 Proposition. If Condition B or C is satisfied (and hence both of them are), then 

EjQ^(u) — > oo if maxUj — minu; — >oo. 

' ' ieV ieV 

3.5.4 Corollary (Solvability alternative). Suppose Condition B or C is satisfied (and hence both 
of them and Condition A are). Then Ej q^(u) (maybe with some variables U; fixed) has a min- 
imizer u^i,^. Either u^in ^ contained in the domain where all triangle inequalities are satisfied, 
then it is unique (up to a constant if no variables are fixed) and corresponds to the solution of 
the corresponding discrete conformal mapping problem, or it lies outside that domain, then the 
corresponding discrete conformal mapping problem does not have a solution. 

Proof of Proposition 3.5.3. Using the (constant) angles a we can rewrite the sum over vertices 
on the right-hand side of equation (3.7) as a sum over triangles: 

iev ijkeT 

Expressing u in terms of A and A, we obtain 

+ i^-KXhi-^ki), 
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Figure 9. The sine theorem lets us express the length-cross-ratios (see Figure 3) in terms angles. 



so 



ijkeT 

where here and in the following ''const." stands for terms that do not depend on u. The estimate 



yields 



f{x,y,z) > n max{x, y,z} 



^T,e,A(") ^ minfaf^.} ^ {max{Xij,Xj^,Xki} - mm{Xij,Xjj^,Xj,i}) + const. 

ijkeT 



Now if ijk e T, then Uj — Uj = Xj^i — Ajj. — Aj-; + ?^jk, so 

max{ui,uj,u;t} -min{Ui,"j,"fc} < max{A;j,Aj;t.-^fci} -min{A;j, A^it-^fci} 

and because the triangulated surface is connected this implies 

£t« j(u) > minfa^lfmaxu; — minu.l + const. 



■ const. . 



D 



3.6 The second variational principle 

The second variational principle has angles as variables. It is based on the two elementary 
observations that, first, the sine theorem lets us express the length-cross-ratios in terms of 
angles, 

sin(a^,)sin(a' ) 



lcr;j 



(3.10) 



sin(a|psin(a[;) 

(see Figure 9), and that, second, if we know the angles in a euclidean triangulation, then we 
can (again using the sine theorem) reconstruct the lengths up to a global scale factor. 
For a triangulated surface T and A e M^, define 



'T,A 



ijksT 
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Figure 10. Left: Triangulation T (solid) and the corresponding graph F (dotted). Right: The 
cycle of edges of F (dashed) that corresponds to an interior edge ij e Ej, and the corresponding 
tangent vector Vj,- e M.^ to Cq. 



where 

V(a,^,r) = JI(a) + JI(^) + JI(r). (3.12) 

3.6.1 Proposition (Rivin [29]). The function V is strictly concave on the domain 

{(a,/3,r)e(K>o)^|a + ^ + r = 7i}. 

So Sjp, is also strictly concave on the domain of positive angle assignments that sum to n in each 
triangle. 



3.6.2 Proposition (Second variational principle I). Let "iS be a discrete conformal class on T 
with representative i = e^l"^ e M^, let 9 e M^, and define the subset C c.M^by 



\ a I 



a > 0, for all ijk e T: a' + ai. + a^. 

JK /CI Ij 

foralli^V: ^ (^jk = ®i} 



n. 



(3.13) 



jk:ijkeT 



Then a & Cq is the angle function of a solution I = e^^^ of Problem 3.1.1 if and only if Sjx{.a) 
is the maximum of the restriction Sj p^L . 

Proof. Consider the graph F that is obtained by choosing one point in each triangle of T and 
connecting it to the vertices of the triangle (see Figure 10, left). The vertex set Vp is in one-to- 
one correspondence with Vj U Tj, and the edge set Ep is in one-to-one correspondence with 
the set of angles Aj. The tangent space to Cq Q R^, which consists of those vectors in W^ that 
sum to in each triangle and around each vertex, is in one-to-one correspondence with the 
space of closed edge chains of F with coefficients in M. 

First, assume that a is a critical point of Sj^lc- Suppose ij e Ej is an interior edge and 
consider the cycle of edges of F shown on the right in Figure 10. The corresponding tangent 
vector to Cc, in M^ is 



e 



dal, 



da]. da 




jk 



da 



(3.14) 



ki 
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Figure 1 1 . Edge-connected chain of triangles from edge ij to edge kl of T. The corresponding 
chain y ofF-edges (dotted) consists of the edges ofF opposite the initial edge ij, the terminal edge 
kl, and the intermediate connecting edges of the triangle chain. 



and 



l^ - ^ + ^ - ^J^T,A(a) = -log 1^-— — — 



+ Ki -^u + ^ik -^ki 



jk 



Provided that a is in fact the system of angles of a discrete metric £ = e^'^, this implies that 



= e^l\ 

i and £ are discretely conformally equivalent. It remains to show^ that a is indeed the system 
of angles of a discrete metric. Construct such a metric as follow^s: Pick one edge ij e Ej and 
choose an arbitrary value for X^j. To define X^^ for any other edge Ik e Ej, connect it to ij by 
an edge-connected sequence of triangles as shovwi in Figure 11, let 7 be the chain of edges of 
the graph F as indicated in the figure. Let w^ be the corresponding vector in M'^ and define 



^ki — ^ij + dSj 



(Wy) + A 



kl 



A;; 



The value of Aj.; obtained in this w^ay is independent of the choice of triangle chain: Another 
triangle chain connecting ij to kl leads to an edge-chain y' such that 7' — 7 is a closed edge- 
chain so that Wyf — w^ ^ R^ is tangent to Cq. Further, I = e^'^ is a discrete metric with 
angles a: If ij and kl belong to the same triangle (that is, if the triangle chain consists of only 
one triangle) then this follow^s from the sine rule. The general case follow^s by induction over 
the length of y. So £ is a solution of Problem 3.1.1. 

The converse implication (solution of Problem 3.1.1 implies critical point) follows from 
the fact that the cycle space of F is spanned by the cycles corresponding to interior edges of T 
as shown in Figure 10 (right) together with the cycles in F corresponding to edge-connected 
triangle sequences as shown in Figure 11 but closed. D 

Solutions to the more general Problem 3.1.3 (with u\y = 0) are also in one-to-one cor- 
respondence with critical points of Sjx- The only difference is that the angle sums are not 
constrained for vertices in Vq- 
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3.6.3 Proposition (Second variational principle II). Let '^ he a discrete conformal class on T 
with representative I = e'^^^ e R^, let V = VqU^i ^^ a partition of V, let Q e M^i, and define 

Cq c r^ by 



Ce 



e = {«' 



a> 0, for all ijk e T: a}., + a^,. + a^. = n, 

]K Kl IJ 

foralli^V,: J] (^jk = ®i}- 

jk:ijkeT 



(3.15) 



Then d ^ Cq is the angle function of a solution of Problem 3.1.3 with fixed u\y = if and only 
ifSj x(ct) is the maximum of the restriction Sj x\r' ■ 

We omit the proof because no essential new ideas are necessary beyond those used in the 
proof of Proposition 3.6.2. 

3.6.4 Remark. It also makes sense to consider critical points of Sj^ under variations of the 
type shovwi in Figure 10 alone, disallow^ing variations corresponding homologically non- trivial 
cycles in F. These correspond to discretely conformally equivalent similarity structures, that is, 
"metrics" which may have global scaling holonomy 

3.7 Mapping to the sphere and to the disk 

If one can solve Problem 3.1.3, one can also find discrete conformal maps from euclidean 
triangulations that are topological spheres to polyhedra with vertices on the unit sphere in 
M.^, and from euclidean triangulations that are topological disks to planar triangulations with 
boundary vertices on the unit circle. The method for mapping to the sphere presented here is 
a slight variation of a method that was previously described [34]. The old version requires an 
input triangulation that is immersed in some M" with straight edges. The variation presented 
here requires only a triangulated surface with discrete metric (but yields the same maps for 
immersed triangulations). 

Mapping to the sphere. Suppose (T, £) is a euclidean triangulation that is topologically a 
sphere. Proceed as follows: 

1. Choose a vertex k and apply a discrete conformal change of metric (2.1) so that afterwards 
all edges incident with k have the same length, say £j-; = 1 for all neighbors i of k. For 
example, let e"'/^ = l^^ if i is a neighbor of k and 1 otherwise. 

2. Let T' be T minus the open star of k. This is topologically a closed disk. 

3. Solve Problem 3.1.3 for T' with prescribed 0; = 2n for interior vertices i and prescribed 
Uj = for boundary vertices. (Suppose a solution exists.) The result is a planar triangula- 
tion. 

4. Map the vertices of this planar triangulation to the unit sphere by stereographic projection. 
Add another vertex (the image of the removed vertex fc) on the sphere at the center of the 
stereographic projection. Build a geometric simplicial complex using these vertices and the 
combinatorics of T. 

3.7.1 Proposition. The result of this procedure is a polyhedron with vertices on the sphere that 
is discretely conformally equivalent to (T, £). (It may not be convex. It is also possible that the 
planar triangulation obtained in step three overlaps with itself In this case the star of k in the 
image polyhedron is not embedded.) 
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Figure 12. Mapping to the disk. 



Proof. After Step 1, the length-cross-ratio for an edge ki incident with k is the quotient of the 
lengths of two consecutive edges mi, ij in the boundary of T'. This is not changed in Step 3 
because u = on the boundary. Further, the length-cross-ratio for an edge ij opposite k as in 
Figure 3 is then the quotient ia/iij- This is also not changed in Step 3 because u = on the 
boundary. Now imagine that before Step 4 you reinsert fc at oo in the plane, which you identify 
with the (extended) complex plane. Then the absolute values of the complex cross-ratios for 
all edges are the same as in (T, £). D 

Mapping to the disk. Suppose (T, £) is a euclidean triangulation that is topologically a closed 
disk. Proceed as follows (see Figure 12): 

1. Choose a boundary vertex k and apply a discrete conformal change of metric (2.1) so that 
afterwards all edges incident with k have the same length, say Ij^j = 1 for all neighbors i of 
k. For example, let e"''^ = £^.^ if i is a neighbor of k and 1 otherwise. 

2. Let T' be T minus the open star of k. Suppose this is topologically still a closed disk. 

3. Solve Problem 3.1.3 for J', with prescribed 9; = 2n for interior vertices of T', 9; = ti 
for boundary vertices of T' that are not neighbors of k in T, and prescribed U; = for the 
neighbors of k in T. (Suppose a solution exists.) The result is a planar triangulation. All 
boundary edges except the neighbors of fc in T are contained in one straight line. 

4. Apply a Mobius transformation to the vertices that maps this straight line to a circle and 
the other vertices inside this circle. Reinsert k at the image point of oo under this Mobius 
transformation. 

3.7.2 Proposition. The result of this procedure is a planar triangulation that is discretely con- 
formally equivalent to (T, £) and has a boundary polygon that is inscribed in a circle. 
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Figure 13. Parallel transport of directions from edge to edge. 

We omit the proof because no new ideas are needed. 

3.7.3 Remark. Note that it is a necessary condition for the solvability of Problem 3.1.3 in 
step 3 that T has no ears (triangles with two edges on the boundary). 

3.8 On the corresponding smooth theory 

A natural question regarding the two variational principles for discrete conformal maps is: 
"What are the corresponding variational principles in the classical smooth theory of conformal 
maps?" In fact, even the question "What exactly are the corresponding smooth mapping prob- 
lems?" deserves a comment. For the second variational principle it is not even obvious how 
the variables — triangle angles — translate to the smooth theory. 

To answer this last question, consider a system a e M'^ of angles in a triangulation T 
that sums to n in each triangle. Even if this angle assignment does not allow a consistent 
assignment of edge lengths, it does allow a sensible definition of parallel transport of directions 
from edge to edge: The direction that makes an angle /3jj with the directed edge ij in triangle 

ijk makes an angle /3jj, = p^j + a^. — n with edge jk (see Figure 13). So angle assignments in 
triangulated surfaces correspond to connections in the direction bundle of a smooth surface. 

Our discussion of the second variational principle will focus on the special case of flat 
connections in the direction bundle (this corresponds to angle sum 2n around interior vertices) 
that also have trivial global holonomy For any such connection there exists a parallel direction 
field and this is unique up to rotation by a constant angle. Conversely, any direction field is 
parallel for a unique flat connection with trivial global holonomy. 

Before we enter our discussion of the smooth conformal mapping problems and variational 
principles, we will first review some background material about conformal metrics on Riemann 
surfaces from the point of view of direction fields. This is of course classic material, only 
the presentation is a bit unusual. Nevertheless, the smooth version of the second variational 
principle seems to be little known, although it is essentially just the theory of harmonic maps 
from a surface to S^. 

Let M be a smooth oriented surface, possibly with boundary, equipped with a Riemannian 
metric g. Let J : TM — > TM be the rotation by 90°, so that the area 2-form is cr = g(,J-, •), and 
let V be the Levi-Civita connection of g. The tensor J is parallel, that is, Vx('^^) — ■^C^x^)- 
The Riemann curvature tensor is 

K(X, Y)Z = VxVyZ - VyW^Z - S7[x,Y]Z, 

where [•, •] is the Lie bracket for vector fields. 

By a direction field on M with respect to g we mean simply a vector field Y with g(Y, Y) = 1, 
that is, a unit vector field with respect to g. Of course the existence of a direction field imposes 
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certain restrictions on the topology of M. In any case, direction fields exist locally, so purely 
local considerations remain valid for arbitrary M. 
The Gauss curvature of g is 

K=-gliRliY,JY)Y,JY), 

for any direction field Y. The curvature 2-form of g is 

n = Ka. 

Define the rotation 1-form py of a direction field Y with respect to g by 

PY(X) = g{_VxY,JY). 

3.8.1 Proposition. For any direction field Y, 

dpy = —D.. 
Proof. The claim follow^s by a straightforw^ard calculation: 

dpy{Y,JY) =Y ■ giVjyYJY) - {JY) ■ giVyYJY) - g{y^y^jy^Y,JY) 

=g(Vy Vji^y, jy) + giVjyY, VyjY) 

=0 

- g(Vji,Vyy, jy) - g(VyY,VjyJY) -g(V[yjy^Y,JY) 

=0 

=g(i?(y,jy)y,jy) = -ic. 

We have used that VyY _L Y and VyJY _L JY for any vector field V, so that, because M is 
two-dimensional, gCVyF, VyyJY) = for any V, W. D 

Now consider a conformal change of metric by equation (1.1), which is also characterized 
by the fact that the 90° -rotation with respect to g is the same J. The Levi-Civita connection V 
of g is related to the Levi-Civita connection V of g by 

V^Z = VxZ + g(X, G)Z + gliZ, G)X - gliX,Z)G, (3.16) 

where G = grad u, that is, du = g(G,-). A direction field Y with respect to g naturally 
determines a direction field 

y = e""y 

with respect to g. Its rotation 1-form is 

Py(X) = gliVxY,JY). 

3.8.2 Proposition. The rotation 1-forms py, pf are related hy 

Py = Py + *du, 
where * denotes the Hodge star operator for g. 
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The Hodge star operator maps a 1-form o) to the 1-form *a) = —co(J-). It also maps a function 
(0-form) / to the 2-form *f = fa and vice versa, **/ = /. Note that on a 2-dimensional 
manifold, the Hodge star operator for 1-forms depends only on the conformal class of the 
metric. 

Proof of Proposition 3.8.2. First, 

V;fy = e-"(-du(x)y + V;fy), 

so 

Using (3.16) one obtains 

Py(X) =py(X) + g(y, G)gliX,JY) - g(X, Y)gliG,JY). 

Finally, since J is skew, J^ = — 1, and (Y,JY) is an orthonormal frame, 

g(y, G)g(X, jy) - g(X, Y)gliG,JY) = 

-gliG,Y)gliY,JX)-gliG,JY)gliJY,JX) 

= -gliG,JX) = -duUX) = *du(X), 

and this completes the proof. D 

As a corollary of Propositions 3.8.1 and 3.8.2, we obtain the equation relating the curvature 
2-forms of g and g, 

n = n-d*du, (3.17) 

and hence Liouville's equation for the curvatures, 

e^''k = K + Au, (3.18) 

where A is the Laplace-Beltrami operator with respect to g, 

A/ = -*d*df. 

(We use the sign convention for the Laplace operator that renders it positive semidefinite.) 

Which problems in the smooth theory are analogous to the discrete conformal mapping 
problems discussed in this paper? There are two fairly obvious candidates: 

3.8.3 Problem. Given (M, g) and the function K, find a conformally equivalent Riemannian 
metric g with curvature K. 

This amounts to solving Liouville's equation (3.18) for u. 

3.8.4 Problem. Given (M, g) and the 2-form Q, find a conformally equivalent Riemannian 
metric g with curvature 2-form Q. 
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This amounts to solving equation (3.17), which is equivalent to Poisson's equation 

Au=f 

with right-hand side f = *Q — K. 

For both problems, one may prescribe either u on the boundary (Dirichlet conditions) or 
*du\j-g]y[ with J *du = — j (Cl — Q) (Neumann conditions). 

Both Liouville's equation (3.18) and the Poisson equation (3.17) are variational, with the 
respective functional 

£^(u)= (iduA*du-(|e2"K:-uiC)cr) (3.19) 

Jm 

for Liouville's equation, where K and K are fixed given functions on M, and 

£g(u)= {^duA*du-u(n-n)) (3.20) 

Jm 

for equation (3.17), where ri and Q are fixed given 2-forms on M. 

We would argue that the discrete conformal mapping problems in Section 3.1 are analo- 
gous to Problem 3.8.4 rather than to Problem 3.8.3, and that the function fj.e.A of the first 
variational principle corresponds to Eg rather than £^. Although we did denote the angle de- 
fect at vertex i by iC; (Section 3.5), prescribing the angle defect is analogous to prescribing the 
curvature 2-form rather than the curvature, the latter being an angle defect per surface area. 
Probably the clearest manifestation of this difference is in the scaling behavior: The angle 
defects at vertices and the curvature 2-form do not change when all lengths are scaled by a 
constant factor. The curvature on the other hand is inversely proportional to the square of this 
factor. Note, however, that for K = 0, Q = there is no difference between Problems 3.8.3 
and 3.8.4, and Ej^ = Eg. 

The remaining question is: What is the smooth version of the second variational principle 
(Section 3.6)? 

Assume for now that M is diffeomorphic to a closed disk and Jl = 0, that is, we are looking 
for a conformally equivalent flat metric. We will discuss this case in detail, because it allows a 
more intuitive treatment using direction fields rather than connections or connection 1 -forms. 
Afterwards, we will indicate how to treat the general case. 

For a direction field Y define 



Jm 



IM 

3.8.5 Proposition (First variation). Let Y^ be a variation of the direction field Y with 

d 



where h e C°°{M). Then 



d 
dt 



dt 



s(n) = 

t=o 



Ft = hJY, 



h d*pY + 
Jm 



h *p . 

dM 
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Proof. This follows from — I .py = dh. D 

•^ at \t=0 ^ 't 

3.8.6 Corollary. A direction field Y is a critical point of S under variations that fix Y on the 
boundary dM if and only if 

d*pY = 0. 

It is also a critical point of S under arbitrary variations if and only if additionally, 

*Py\t3m = ^- ^3.21) 

Loosely speaking, the following proposition says that straightest direction fields with re- 
spect to g are parallel direction fields with respect to a conformally equivalent flat metric g 
with trivial global holonomy. 

3.8.7 Proposition (Smooth version of the second variational principle). Suppose the direction 
field Y is a critical point of S under variations that fix Y on the boundary. Define the function u 
up to an additive constant by 

du = *Py- 

(This is possible since *Py is closed by the above corollary and we assumed that M was diffeomor- 
phic to a disk.) Define g by the conformal change of metric (1.1). Then: 

(i) The direction field Y = e~"Y is parallel with respect to g, so g is flat. 

(ii) The geodesic curvature of the boundary dM with respect to g is 

k = K - Py{T), 

where k is the geodesic curvature with respect to g and T is the positively oriented unit tangent 
vector field to dM. 

(Hi) IfYis also a critical point of S under arbitrary variations, then u\gf^ is constant. 

Proof. Since *du = — py, the rotation form of pf vanishes by Proposition 3.8.2. This im- 
plies (i). The geodesic boundary curvatures are k = pj(T) and k = pff. (Locally extend 
the direction field T inwards from the boundary.) Again by Proposition 3.8.2, this implies (ii). 
Finally, (iii) follows immediately from (3.21). D 

In the general case, M is not restricted to be diffeomorphic to closed disk and one is looking 
for a conformally equivalent metric g with prescribed curvature 2-form ti. To treat this case 
variationally consider the functional 



S{p) = ] 



2 

Jm 



p A*p 



on the affine space of 1-forms p satisfying dp = (i — Q. We leave the details to the reader, not 
because they are tedious but because they are interesting. The critical points correspond to 
conformally equivalent similarity structures, that is, "metrics" which may have global scaling 
holonomy. (Compare the remark at the end of Section 3.6.) 
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4 Ideal hjrperbolic polyhedra 

4.1 Hj^erbolic structure on a euclidean triangulation 

This section deals with the inverse of a construction of Penner [27] [7], which equips a hy- 
perbolic manifold with cusps with a piecewise euclidean metric. Here, we construct a natural 
hyperbolic metric with cusps on any euclidean triangulation. 

Consider a euclidean triangle with its circumcircle. If we interpret the interior of the cir- 
cumcircle as a hyperbolic plane in the Klein model, then the euclidean triangle becomes an 
ideal hyperbolic triangle, that is, a hyperbolic triangle with vertices at infinity. This construc- 
tion equips any euclidean triangle (minus its vertices) with a hyperbolic metric. If it is per- 
formed on all triangles of a euclidean triangulation (T, £), then the hyperbolic metrics induced 
on the individual triangles fit together so T \ V is equipped with a hyperbolic metric with cusps 
at the vertices. Thus, T becomes an ideal triangulation of a hyperbolic surface with cusps. 

4.1.1 Theorem. Two euclidean triangulations (T, £) and (T, £) with the same combinatorics are 
discretely conformally equivalent if and only if the hyperbolic metrics with cusps induced by the 
circumcircles are isometric. Discrete conformal maps are isometries with respect to the induced 
hyperbolic metrics. 

Proof. This follows immediately from Theorem 2.7.2 (Section 2.7), because the projective 
circumcircle preserving maps between triangles are precisely the hyperbolic isometries. D 

4.1.2 Remark. Each discrete conformal structure on T corresponds therefore to a point in 
the classical Teichmiiller space 5^ „ of a punctured surface. This explains the dimensional 
agreement observed in Section 2.5. 

Theorem 4.1.1 also suggests a way to extend the concepts of discrete conformal equiva- 
lence and discrete conformal maps to triangulations which are not combinatorially equivalent: 

4.1.3 Definition. Two euclidean triangulations (T, £) and (T, £), which need not be combina- 
torially equivalent, are discretely conformally equivalent if they are isometric with respect to 
the induced hyperbolic metrics with cusps. The corresponding isometries are called discrete 
conformal maps. 

This idea will be pursued further in the sequel [33]. Here, we will stick with Defini- 
tion 2.1.1. 

4.2 Decorated ideal triangles and tetrahedra 

In this section we review some basic facts about ideal triangles and tetrahedra that will be 
needed in subsequent sections. 

All ideal hyperbolic triangles are congruent with respect to the group of hyperbolic isome- 
tries. A decorated ideal triangle is an ideal hyperbolic triangle together with a choice of horo- 
cycles, one centered at each vertex (see Figure 14). We denote by A;^ the signed distance 
between the horocycles at vertices i and j as measured along the edge ij and taken negatively 
if the horocycles intersect. Any triple of real numbers (Ajj, A^j^, Aj.;) e M.^ corresponds to one 
and only one choice of horocycles. Figure 14 shows also the lines of symmetry of the ideal 
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Figure 14. Decorated ideal triangle in the Poincare disk model. 



triangle. (They are its heights as well.) We denote the signed distances from their base points 
to the horocycles by p^- as shown. Clearly, 



^ij=pik+Pw 



^jk=pii+ptp 



hi=Pij+p]k' 



so 



Pij — 2^~^ij + ^jk + ^ki)> 
-(A;- - ^jk + ^ki)' 



Jk 



(4.1) 



Pki ~ S^'^O' + ^jk ~ '^fci)- 



4.2.1 Lemma (Penner [27]). The length c^- of the part of the horocycle centered at v^. that is the 
ideal triangle as shown in Figure 14 is 

Proof See Figure 15, which shows the ideal triangle in the half-plane model. Recall that in 
the half-plane model, the hyperbolic plane is represented by {z e C| Imz > 0} with metric 



ds = -^\dz\. 

Imz ' ' 



D 



4.2.2 Remark. Together with Proposition 4.3.2 of the next section, this provides a geometric 
interpretation for the auxiliary parameters c'.^ introduced in Section 2.4. 

All ideal tetrahedra are not isometric. There is a complex 1 -parameter family of them, the 
parameter being the complex cross-ratio of the vertices in infinite boundary of hyperbolic 3- 
space. A decorated ideal tetrahedron is an ideal hyperbolic tetrahedron together with a choice of 
horospheres centered at the vertices. Figure 16 shows a decorated ideal tetrahedron, truncated 
at its horospheres, in the half-space model. Again, we denote the signed distances between 
the horospheres by Aj^. 

The intrinsic geometry of a horosphere in hyperbolic space is euclidean. So the intersection 
of the tetrahedron with the horosphere centered at, say, V; is a euclidean triangle with side 
lengths c| ., c'^^, c[. determined by Lemma 4.2.1. One easily deduces the following lemma. 
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Figure 15. Decorated ideal triangle in the half-plane model. 




Figure 1 6. Decorated ideal tetrahedron in the half-space model. 
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Figure 1 7. Penner coordinates (left) and shear coordinates (right). 



4.2.3 Lemma. Six real numbers X^p Xj^, Xj^^, Xq, X^, Aj.; are the signed distances between 
horospheres of a decorated ideal tetrahedron as shown in Figure 1 6 (which is then unique) if and 
only if c\; c'.^,, c[. determined by Lemma 4.2.1 satisfy the triangle inequalities. 

So the six parameters X determine the congruence class of the ideal tetrahedron (2 real 
parameters) and the choice of horospheres (4 parameters). 

Note that the angles of the euclidean triangles in which the tetrahedron intersects the four 
horospheres are the dihedral angles of the tetrahedron. This implies that the dihedral angles 
sum to 71 at each vertex, and further, that the dihedral angles at opposite edges are equal. 
The space of ideal tetrahedra is therefore parametrized by three dihedral angles a^j = a^.i, 
("■jk = ct;;, a^i = aji satisfying a^ + a^; + a^i 



71. 



4.3 Penner coordinates and shear coordinates 

In Section 4.1, we equipped a euclidean triangulation (T, £) with a hyperbolic cusp metric that 
turns it into an ideal hyperbolic triangulation. In this section, we will identify the logarith- 
mic edge lengths X (see equation (2.2)) with the Penner coordinates [27] and the logarithmic 
length-cross-ratios loglcr (see equation (2.4)) with the shear coordinates [10] [37] for this 
ideal triangulation. (The handbook [25] is a good reference for the pertinent aspects of Teich- 
miiller theory.) 

Warning. Our notation differs from Penner's in a potentially confusing way. His "lambda- 
lengths" are -/2e'^'^ = -/2£ in our notation. Our As are the signed hyperbolic distances be- 
tween horocycles. 

Since the sides of an ideal hyperbolic triangle are complete geodesies, there is a one- 
parameter family of ways to glue two sides together. Penner coordinates and shear coordinates 
can be seen as two ways to describe how ideal triangles are glued together along their edges 
to form a hyperbolic surface with cusps. 

Suppose T is a triangulated surface and X e M^. For each triangle ijk e T, take the 
decorated ideal triangle with horocycle distances X^j, Xj^, X^i, and glue them so that the 
horocycles fit together (see Figure 17, left). The result is a hyperbolic surface with cusps at 
the vertices, together with a particular choice of horocycles centered at the cusps. In this way, 
the Penner coordinates X parametrize the decorated Teichmilller space, that is, the space of 
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Figure 18. Shear and length-cross-ratio (see proof of Proposition 4.3.2). 

hyperbolic cusp metrics on a punctured surface (modulo isotopy) with horocycles centered at 
the cusps. 

The shear coordinates represent another way to prescribe how ideal triangles are glued, 
for which no choice of horospheres is necessary. The shear coordinate Z on an interior edge of 
an ideal triangulation is the signed distance of the base points of the heights from the opposite 
vertices (see Figure 17, right). The following relation between Penner coordinates and shear 
coordinates is well known. 

4.3.1 Lemma, if A e M^ are the Penner coordinates for an ideal triangulation with a particular 
choice of horocycles, then the shear coordinates Z e M^'"f are 

where k and I are the vertices opposite edge ij as in Figure 1 7. 
Proof. The claim follows from Z;^ = p^. — pl^ and equations (4.1). 



n 

lEint fgj- 

(4.2) 

M^ are 



4.3.2 Proposition. Let (T,i) be a euclidean triangulation. The shear coordinates Z e 
the corresponding ideal triangulation (see Section 4. 1) are 

Zij = loglcVij 

(see equation (2.4)j. Thus, for a suitable choice of horocycles, the Penner coordinates A i 
given by equation (2.2). 

Proof. Consider an interior edge ij e E between triangles ijk and jil. Without loss of gen- 
erality, we may assume that the triangles have a common circumcircle. For otherwise we 
can change £ discretely conformally so that this holds, and this changes neither Icr^^ nor the 
hyperbolic cusp metric on T. We may further assume that ij is a diameter of the common cir- 
cumcircle. For otherwise we may apply a projective transformation that maps the circle onto 
itself so that this holds. This is an isometry with respect to the hyperbolic metric of the Klein 
model, and it is a discrete conformal map of the quadrilateral formed by the two triangles. We 
arrive at the situation shown in Figure 18 in the Klein model. The hyperbolic heights are also 
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the euclidean heights, and in the hyperbolic metric of the Klein model, the distance between 
their base points a and b is 

\avj\\bvi\ 

^ 2 \avi\\bvj\ 

where |xy| denotes the euclidean distance between x and y. Since by elementary euclidean 
geometry 

— and — 



this implies equation (4.2). Now Lemma 4.3.1 implies the statement about Penner coordinates. 

D 

4.4 Ideal hj^erbolic polyhedra with prescribed intrinsic metric 

The discrete conformal mapping problems described in Section 3.1 are equivalent to problems 
involving the polyhedral realization of surfaces with hyperbolic cusp metrics, like the follow- 
ing. 

4.4.1 Problem. Given an ideal triangulation T of a punctured sphere equipped with a hy- 
perbolic metric with cusps, find an isometric embedding of T as ideal hyperbolic polyhedron 
in H^. The polyhedron is not required to be convex, but it is required that the edges of the 
polyhedron are edges of T. 

4.4.2 Theorem. For any vertex I of J, Problem 4.4. 1 has at most one solution that is star-shaped 
with respect to I. 

Problem 4.4.1 is equivalent to (a special case of) Problem 3.1.3, so Theorem 4.4.2 follows 
from Theorem 3.1.5. Indeed, to solve Problem 4.4.1 one may proceed as follows. Let A e M'' 
be the Penner coordinates for the ideal triangulation T, and let I = e^'^. Choose a vertex I 
of T and let the triangulation T' be T minus the open star of I. Solve Problem 3.1.3 for J', 
prescribing 0j = 271 if i is an interior vertex and U; = — A;; if i is a boundary vertices. Suppose 
a solution u e M^ exists. This leads to a flat triangulation (T',£). Suppose further that (T',l) 
does not overlap with itself when developed in the plane. For each triangle ijk of J', construct 
the decorated ideal tetrahedron (see Section 4.2) with horosphere-distances ?i^j, Xj^, X^^, and 
—Uj^, —Uj, —U]^ as shown in Figure 19. They exist by Lemma 4.2.3, because by Lemma 4.2.1, 
the intersection of the ideal tetrahedron with the horosphere centered at the vertex opposite 
ijk is the euclidean triangle with side lengths l^j, l^j,, 1^^; as shovwi in the figure. Hence, all 
these ideal tetrahedra fit together to form a solution of Problem 4.4.1 that is star-shaped with 
respect to I. 

Conversely, any solution of Problem 4.4.1 that is star-shaped with respect to I yields a 
solution without self-overlap of the corresponding instance of Problem 3.1.3. 

4.4.3 Remark. Note the similarity with the procedure for mapping to a sphere described in 
Section 3.7. 

Numerous variations of Problem 4.4. 1 can be treated in similar fashion. We mention only 
the following. 
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Figure 19. Discretely conformally flattening a euclidean triangulation is equivalent to construct- 
ing an ideal polyhedron with prescribed hyperbolic cusp metric. 

AAA Problem. Given an ideal triangulation T of a punctured torus equipped with a hyper- 
bolic metric with cusps, find an isometric embedding of the universal cover of T as an ideal 
polyhedron in H^ that is symmetric with respect to an action of the fundamental group of T 
by parabolic isometries. The polyhedron is not required to be convex, but it is required that 
the polyhedron is star-shaped with respect to the ideal fixed point of the parabolic isometries 
and that the edges of the polyhedron are edges of T . 

4.4.5 Theorem. If Problem 4.4.4 has a solution, it is unique. 



4.5 The variational principles and hjrperbolic volume 

The connection with hyperbolic polyhedra elucidates the nature and origin of the variational 
principles for discrete conformal maps. In this section, we will indicate how to derive the first 
variational principle for discrete conformal maps (Proposition 3.3.3) from Milnor's equation 
for the volume of an ideal tetrahedron and Schlafli's formula. 

Milnor [22] [23] showed that the volume of an ideal tetrahedron with dihedral angles a, 
13, Y is V(a, /3,y) as defined by equation (3.12). Schlafli's differential volume formula (more 
precisely, Milnor's generalization which allows for ideal vertices [24]) says that its derivative 
is 



1 V^ 
dV = --2_^Xijdaij, 



(4.3) 



where the sum is taken over the six edges ij, X^j is the signed distance between horospheres 
centered at the vertices i and j, and a^j is the interior dihedral angle. (The choice of horo- 
spheres does not matter because the angle sum at a vertex is constant.) Define 



V(Ai2, .^23> ■^31> -^14' -^24' -^34) ~ 2 Z_l "O'^U "*" ^("l4' Ct24' ^'34)' 



(4.4) 
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where the dihedral angles ai2 = 0:34, ot23 = ^1^4, a^i = 024 of the decorated tetrahedron are 
considered as functions of the Aj^. (They are the angles in a euclidean triangle with side lengths 
g(Ai2-Ai4-A24)/2^ g(A23-A24-A34)/2^ gCAjj - A34- Ai4)/2^ ggg Section 4.2.) Then, by equation (4.3), 

dV = -Y,cc,jdXij. (4.5) 

This implies Proposition 3.3.2 on the partial derivatives of Ejq^, and therefore Proposi- 
tion 3.3.3 (the first variational principle), because using (4.4) (and (2.3)) we can rewrite 
equation (3.3) as 

ET,0,xiu) = ^ 2V{Xij, Ajfc, Afc;, -Ui, -uj, -Ufc) 

ijkeT 

Sr . Sr ^4.6) 

ijeE iev 



where 



(See also Figure 19.) 



(71, if i j is an interior edge, 
(4.7) 
^, if ij is a boundary edge. 



5 H5q)erbolic triangulations 

5.1 Discrete conformal equivalence of hjrperbolic triangulations 

In Section 4.5 we derived the first variational principle for discrete conformal maps from Mil- 
nor's equation for the volume of an ideal tetrahedron and Schlafli's formula. A straightforward 
modification of this derivation leads to a companion theory of discrete conformality for hy- 
perbolic triangulations. This makes it possible, for example, to construct discretely conformal 
uniformizations of higher genus surfaces as shown in Figure 20. 

We will present the basic theory in this section, and show how to derive it by hyperbolic 
volume considerations in the next. For the discretely conformal uniformization of a euclidean 
triangulation over the hyperbolic plane, as shown in Figure 20, we need to define what it 
means for a euclidean triangulation to be discretely conformally equivalent to a hyperbolic 
triangulation. This is the topic of Section 5.3. 

Suppose T is a surface triangulation and i e M^q is a discrete metric, that is, a real valued 
function on the set of edges that satisfies all triangle inequalities. Then there is not only a 
euclidean triangulation (T, £) with these edge lengths. One can equally construct hyperbolic 
triangles ijk with hyperbolic side lengths i^j, Ij^., £j,; and glue them together. The result is a 
hyperbolic triangulation which we denote by {T ,l\. 

5.1.1 Definition. Two combinatorially equivalent hyperbolic triangulations, (T, £)/, and (T, i\, 
are discretely conformally equivalent if the discrete metrics I and I are related by 

£ i 

sinh^ = e2("-+"i)sinh^ (5.1) 

2 2 

for some function u : V — > M. 
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Figure 20. Discretely conformal uniformization of a genus-two surface. 

Thus, in the hyperboHc version of the theory, equation (2.2) is replaced by 

i 
A = 21ogsinh-, (5.2) 

so that in terms of A and A, the relation (5.1) is again (2.3). The role of £j g ;^^(u) is played by 
the function 

^T,e,A(") = Xi 2^h(^i;. ^;fc' ^fci' ~"i' ~^P ~^k, ) + Xi ®i"i ' ^^'^^ 

ijkeT iev 

where 

2Vj,(Ai2, A23, A31, Ai, A2, A3) = 

aiAi + a2A2 + a3A3 + ai2Ai2 + «23'^23 + «31'^31 

+ Jl(ai) + JI(a2) + JI(a3) + JI(ai2) + JI(a23) + JlCagi) 

+ Jl(i(7r-ai-a2-a3)), (5.4) 

and ai, a2, a^ are the angles in a hyperbolic triangle with side lengths 

£23 = 2arsinh (65^^23-^2-^3)^^ 

£31 = 2arsinh(ei^^3i-A3-Ai)^^ (5.5) 

£12 = 2arsinh (ell^^i^-^i-^^))^ 
and 

«23^|(^ + «l-«2-«3). 

«31 = |("-«l + «2-«3). (5.6) 



«12 = K"~"l~"2 + a3)- 
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Thus, Vft is defined only on the domain where £i2> ^23' ^3i satisfy the triangle inequalities. 
However, exactly as in the case of Ejqj^(u), we can extend the domain of definition of 
£^ Q ^(u) to the whole of M^: 

5.1.2 Proposition. Extend the domain of definition of Ej ^ ^(u) to M^ by declaring the angles in 
"broken" triangles to be 0, 0, n, respectively. The resulting function is continuously differentiable 

onR^. 

Proof See Section 5.2. D 

5.1.3 Remark. To compute the angles a, p, y in a hyperbolic triangle with side lengths a, b, 
c, one can use, for example, the hyperbolic cosine rule or the hyperbolic half-angle formula 



r a A _ / sinh ((a - b + c)/2) sinh ((a + b - c)/2) 
^^" y2j~\j sinh ((-a + b + c)/2) sinh ((a + b + c)/2) ' 

5.1.4 Proposition. Let I e M^, let X be defined by equation (5.2), and suppose u e M^ is in 
the domain where I defined by equation (5.1) satisfies all triangle inequalities. Then the partial 
derivative of Ej ^ ^ vifith respect to U; is 

^"' jk:ijkeT 

where a are the angles in the hyperbolic triangulation (T, l)j^, and the sum is taken over all angles 
around vertex i. 

Proof See Section 5.2. D 

5.1.5 Proposition. The function Ej^^i^u) is convex on M^ and locally strictly convex in the 
domain where I defined by equation (5.1) satisfies all triangle inequalities. 

Proof. See Section 5.2. D 

Consider the discrete conformal mapping problems for hyperbolic triangulations that are 
analogous to those for euclidean triangulations described in Section 3.1. Propositions 5.1.4 
and 5.1.7 imply the following hyperbolic version of Theorem 3.1.5. 

5.1.6 Theorem. If the discrete mapping problems for hyperbolic triangulations have a solution, 
it is unique and can be found by minimizing Ej „ ;^(")- 

In contrast to the euclidean case, there is never any ambiguity of scale. The following relatively 

simple explicit formula for the second derivative facilitates the numerical minimization of 

E^ 
T,ex 

5.1.7 Proposition. The second derivative of Ej ^ ^atu is 



d'^E^ 



. . ,, dwdUi 

i,jeV ' J 



]_ Y,w,j{u)[{du,-dUjf + tanh\^-f){du, + du^f'), (5.7) 



ijGE 
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Figure 21. The polyhedral building block (shown in the Poincare ball model) used to derive the 
theory of discrete confomial equivalence of hyperbolic triangulations. 

with i defined by equation (5.1) and 

Wij(u) = - (cot (i(7i - a), - al + af.)) + cot (i(7r - 4 - a\. + aj,))) 

for interior edges ij with opposite vertices k and I if I satisfies the triangle equalities for ijk and 
jil (so that the corresponding angles a are positive and smaller than n). If ij is a boundary 
edge, there is only one cotangent term. For "broken" triangles, replace the three corresponding 
cotangent terms with 0. 

We omit the proof, which consists of a lengthy but elementary calculation. 



5.2 Derivation by volume considerations 

The theory of discrete conformal equivalence for hyperbolic triangulations is based on volume 
considerations for the type of polyhedron shown in Figure 21. From the vertices pi, p2, Ps of 
a triangle in hyperbolic 3 -space, three rays run orthogonally to the plane of the triangle until 
they intersect the infinite boundary in the ideal points v^, V2, V3. The convex hull of these six 
points is a prism with three ideal vertices and right dihedral angles at the base triangle PiP2P3- 
Let the dihedral angles at the three rays from pi, p2, Ps be a^, a2, a^- Since the dihedral 
angles sum to n at the ideal vertices, the dihedral angles ai2, 023, a^i at edges v^V2, V2V3, 
V3V1 satisfy equations (5.6). Let i^j be the lengths of the finite edges, and let Aj and Aj^ be the 
lengths of the infinite edges truncated at some horospheres centered at the ideal vertices V;, as 
shovwi in Figure 21. 

5.2.1 Lemma. The (truncated) edge lengths of the prism shown in Figure 21 are related by 
equations (5.5). 

Proof. We consider the case when Aj = A2 = A3 = (that is, when the truncating horospheres 
touch the base plane in pj, p2, P3), from which the general case follows easily. Figure 22 shows 
one of the side quadrilaterals of the prism in the half-plane model. We will show that 



1^ 
A = 21ogsmh(-J, 



(5.8) 
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Figure 22. Proof of Lemma 5.2.1. 



which proves this special case. We have 



-e'-l 



2e* 



01 = 2arccot(^— — j, 02 = 2arccot (^^ j, 

The equation for the "angle of parallelity" (see Figure 23) implies that 

X = logcot(0i/2) - logcot(02/2), 
and hence equation (5.8). 



n 




Figure 23. The angle (p and side length a in a right angled hyperbolic triangle with an ideal 
vertex satisfy the equation a = logcot(</)/2). 
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The volume of the polyhedron shown in Figure 21 is 

If 

V;,(ai, a2, ag) = 3 (JI(«i) + JI(«2) + JlCag) + JlCa^s) + JI(a23) + JI(«3i) 

+ Jl(|(7r-ai-a2-a3))). (5.9) 

This was shown by Leibon [18], who also showed that the volume function V), is strictly con- 
cave on its domain of definition, 

{(ai, a2, as) e M'' a^ > 0, a2> 0, 03 > 0, a^ + a2 + as < ti}. 

By Schlafli's formula, 

^\ = -| (^1 ^^l + ^2 d«2 + ^3 d«3 + ^12 d«12 + ^23 d«23 + ^31 d«3l) • 

(The choice of horospheres does not matter because the angle sum at the ideal vertices is 
constant. Also note that the lengths i^j of the finite edges do not appear in the equation 
because their dihedral angles are constant.) 
The function V;, defined by equation (5.4) is 

^h('^12> ■^23> -^31> -^1. ■^2. -^3) = 

|(aiAi + a2A2 + agAg + 0^2^12 + a23A23 + aai A31) + Vft(ai, a2, ag), 

so that 

d Vft = I (ai dXi + a^ d A2 + ag d A3 + ai2 d A12 + a23 d A23 + aj,x "i^si) • 

From this one obtains Proposition 5.1.4 on the partial derivatives of £^ ^ ^. By extending V^ 
using essentially the same argument as in the proof of Proposition 3.3.5, one obtains Propo- 
sition 5.1.2 on the extension of £^ „ ^. To prove the convexity of £^ ^ ^^^ (Proposition 5.1.5), 
note that the function 

^h(-^12. -^23. -^31. ■^1. -^2> -^3) ~ f (-^12 + ^^23 + -^3l) 

really only depends on the three parameters 

9V, 



^1 — 4( + '^12-'^23 + '^31 -2Ai) — , 

^2 = 4 ( + ^^12 + ^^23 - ^^31 - 2A2) = 



5y, 



5a2' 

^3 = 4 ( - ■^12 + -^23 + ■^31 - 2A3) = 



dcLz 



As function of these parameters, it is minus the Legendre transform of the strictly concave 
function V^^ : 

%-\ (^12 + ^23 + ^31) = -«1^1 - «2^2 - «3^3 + Vft. 

Hence, VJ, is a locally strictly convex function of x^, x-^, X3, and hence also of Aj, A2, A3, if A12, 
A23, A31 are considered constant. The C^ extension of V^ is linear outside the domain where 
the triangle inequalities are satisfied, hence still convex. 
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5.2.2 Remark. In the same way, one can derive a theory of discrete conformal equivalence for 
spherical triangulations. In this case, the polyhedral building block is a tetrahedron with one 
finite and three ideal vertices. However, the functions involved in the corresponding variational 
principles are not convex. So in this case, the variational principles do not immediately lead 
to a uniqueness theorem, nor to a computational method for discrete conformal maps. 

5.3 When are a euclidean and a hjrperbolic triangulation discretely conformally 
equivalent? 

We propose the following definition: 

5.3.1 Definition. A euclidean triangulation (T, £) and a hyperbolic triangulation (T, £);, are 
discretely conformally equivalent if i and I are related by 

I 
sinh^ = e2^"-+"-'^£;j- (5.10) 

for some function u e M^. 

This is based on the following interpretation of equations (5.2) and (2.3). Consider the 
hyperboloid model of the hyperbolic plane, 

H^ = {x e M^ I (x, x) = -1, X3 > 0}, 
where (•, •) denotes the indefinite scalar product 

(x, y ) = xiji + X2y2 - ^3^3, 
and the hyperbolic distance dj^(x,y) between two points x,y e H^ satisfies 

cosh4(x,y) = -(x,y). 

This implies 

lk-jllh = 2sinhQ4(x,j)J, 

where ||v||ft= y{v,v). 

To every triangle hyperbolic triangle in H^ with sides of length £i2) ^23. ^3i) there corre- 
sponds a secant triangle whose sides are the straight line segments connecting the vertices in 
M'^. Their lengths, as measured by || • ||, are therefore 2sinh(£jj/2). 

5.3.2 Remark. Let Pi,P2,P3 e H^ c M"^ be the vertices of a hyperbolic triangle. Then the 
following statements are equivalent: 

(i) The restriction of the semi-Riemannian metric (•, •) on M.^ to the affine plane in spanned 
by pi, P2, P3 is positive definite and therefore induces a euclidean metric on that plane, 
(ii) The side lengths of the secant triangle satisfy the triangle inequalities, 
(iii) The circumcircle of the hyperbolic triangle is a proper circle. (The circumcircle of a 
hyperbolic triangle is either a proper circle or a horocycle or a curve of constant distance 
from a geodesic.) 
The analogous statement for ideal triangles is Penner's Lemma 2.2 [27]. 

Now let (T, £)ft and (T, £);, be a two hyperbolic triangulations. The edge lengths of the 
secant triangles of (T, £)^ are 2e^^^ with A defined by equation (5.2), and similarly for (T, £);,. 
Now (T,l)h and iT,l)h ^re discretely conformally equivalent if and only if £' = e'^^^ and l' = 
e '^ are related by equation (2.1), that is, related like discrete metrics of discretely conformally 
equivalent euclidean triangulations. 
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Figure 24. Circumcirde intersection angles ^ij = a^. + aL For a boundary edge ij e Ei,^y, define 



^n 



a; 



6 Circle patterns, circular meshes, and circle domains 

6.1 Two variational principles for circle patterns 

While the discrete conformal mapping problems essentially ask for ideal hyperbolic polyhedra 
with prescribed metric, the circle pattern problem below asks for an ideal polyhedron with 
prescribed dihedral angles. Rivin's variational principle for this type of problem is very similar 
to our second variational principle for discrete conformal maps. The function is essentially 
the same, only the constraints placed on the angle assignments are stronger. Also, the first 
variational principle needs only a slight modification to become a variational principle for 
circle patterns. 

6.1.1 Problem (Circle pattern problem). Given a surface triangulation T and a function <l> e 
(0, Ti]^, find a discrete metric I so that the euclidean triangulation (T, £) has circumcirde 
intersection angles <l>jj as shown in Figure 24. 

6.1.2 Proposition (Rivin [29]). Angles d&R^ are the angles of a euclidean triangulation (T, £) 
that solves Problem 6.1.1 if and only ifSjx(Si) is the maximum o/Sj ;^ on the set all a e M^ that 
satisfy 

(i) a>0, 
(ii) a^- + a'.j. + a^. = 7i for all triangles ijk e T, 

(Hi) a^. + a\j^ = ^ij for all interior edges ij e E^„[, 

(iv) a^. = ^ij for all boundary edges ij e E^^^. 

Due to condition (Hi), the choice of the parameter A e M^ of Sj j^ does not matter because 

ijeE 

So in connection with circle patterns, it makes sense to consider only 

Now consider the first variational principle for discrete conformal maps. For <!> e M^, 
e e M^ define 



<§•- 



T,$,e • 



S- 



T,4.,e(^'")= Xi '^^^^ip^jk'hu-Ui,-Uj,-u^)-^^ijXij + ^eiUi. (6.1) 



ijkeT 



ijsE 



isV 
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Compare equation (4.6). If ^ is defined by equation (4.7), then 

So if we fix A and vary u, then we obtain the first variational principle for discrete conformal 
maps. If, on the other hand, we fix u and vary A, then we obtain a variational principle for 
circle patterns (see Proposition 6.1.3 below). Interpret the circumcircles and triangle sides 
as hyperbolic planes in the half-plane model. Then, using equation (4.5) and the fact that 
opposite dihedral angles in an ideal tetrahedron are equal, one gets for an interior edge ij e 
Eint that 

where a are the angles in the euclidean triangulation (T, £) with l determined by equa- 
tions (2.1) and (2.2). (Here we assume that l satisfies the triangle inequalities. Otherwise 
the angles are or ti as stipulated in Proposition 3.4.4.) In the same way, one gets for a 
boundary edge ij 

3 



This implies the following variational principle for circle patterns. (Note that (fj^ ©( A, u) is 
also convex if we fix u and consider A as variables.) 

6.1.3 Proposition. The function i e M^q defined in terms ofX and u by equations (2.1) and (2.2) 
is a solution of Problem 6.1.1 if and only if i satisfies all triangle inequalities and <g'T,e,e(A, u) is 
the minimum of the function A >-> Sjqq^X, u). That is, uis arbitrary but constant. (Without loss 
of generality one could fix u = 0.) 

6.2 Discrete conformal equivalence for circular polyhedral surfaces 

In this section, we generalize the notion of discrete conformal equivalence from surfaces com- 
posed of triangles to surfaces composed of polygons inscribed in circles. The variational prin- 
ciple described below is like a mixture of the first variational principle for discrete conformal 
maps of Proposition 3.3.3 and the variational principle for circle patterns of Proposition 6.1.3. 
An (abstract) polyhedral surface is a surface that is a CW complex. A euclidean polyhedral 
surface is a polyhedral surface obtained by gluing euclidean polygons edge-to-edge. If all of 
the polygons have a circumscribed circle, we speak of a (euclidean) circular polyhedral sur- 
face. A circular polyhedral surface is determined by the polyhedral surface P and the function 
I e (M>o)^'' that assigns to each edge its length. Conversely, a function i e (M>o)^'' defines a 
circular polyhedral surface if and only if it satisfies the "polygonal inequalities": In each poly- 
gon, the length of any edge is smaller then the sum of lengths of the other edges. If i satisfies 
these conditions, we denote the resulting circular polyhedral surface by (P,£). 

6.2.1 Definition. Two circular polyhedral surfaces, (P,£) and (P,£), are discretely conformally 
equivalent if £ and £ are related by equation (2.1) for some function u e M^p. 

To solve the discrete conformal mapping problems for circular polyhedral surfaces that are 
analogous to those described in Section 3.1, proceed as follows: First triangulate the non- 
triangular faces of the given circular polyhedral surface (P,£) to obtain a euclidean triangula- 
tion (T, I) (where I : £j — > M>o. (-Iep — O- Then define <I> by equations (4.7) and minimize 
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Sj,pq[X, u), where X^ = 21og£jj is held fixed if ij e Ep and considered a variable if Ej\Ep, 
and the U; are variables or fixed depending on the mapping problem, as in the case of trian- 
gulations. If £ determined by equations (2.1) and (2.2) for the minimizing (A, u) satisfies the 
triangle inequalities, it is a solution of the mapping problem. 

Note that the values l^j for edges ij e £j \ £p do not enter because the corresponding X^j 
are variables. 

6.3 Discrete circle domains 

A domain in the Riemann sphere C is called a circle domain if every boundary component is 
either a point or a circle. Koebe conjectured that every domain in C is conformally equivalent 
to a circle domain. For a simply connected domain, this is just the Riemann mapping theorem. 
Koebe himself proved the conjecture for finitely connected domains, and after various gener- 
alizations by several other people. He and Schramm gave a proof for domains w^ith at most 
countably many boundary components [13]. (Their proof is based on circle packings.) 

The method for mapping to the sphere described in Section 3.7 w^orks (mutatis mutandis) 
also for the circular polyhedral surfaces discussed in the previous section. This allow^s us to 
map euclidean triangulations to "discrete circle domains", that is, domains in the plane that 
are bounded by circular polygons. 

Suppose (T, £) is a euclidean triangulation that is topologically a disc with holes. To map 
(T, £) to a discrete circle domain, simply fill the holes by attaching a face to each boundary 
polygon and map the resulting circular polyhedral surface to the sphere. 

Note that for a topological disk with holes, we recover in a different guise the procedure 
for mapping to a disk that was described in Section 3.7. 
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